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Abstract 

The problem of cosmological graviton creation for homogeneous and isotropic 
universes with elliptical (e = +1) and hyperbolical (e = —1) geometries is ad- 
dressed. The gravitational wave equation is established for a self-gravitating 
fluid satisfying the barotropic equation of state p = (7 — 1)/?, which is the source 
of the Einstein's equations plus a cosmological A-term. The time dependent 
part of this equation is exactly solved in terms of hypergeometric functions for 
any value of 7 and spatial curvature e. An expression representing an adiabatic 
vacuum state is then obtained in terms of associated Legendre functions when- 
ever 7 7^ § (2n-i) ' wnere n * s an integer. This includes most cases of physical 
interest such as 7 = 0, 4/3,1. The mechanism of graviton creation is reviewed 
and the Bogoliubov coefficients related to transitions between arbitrary cosmic 
eras are also explicitly evaluated. 
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1 Introduction 



After the pioneer works of Grishchuk |T|], the formation of a gravitational wave (GW) 
background in the course of the cosmological expansion has been extensively studied 
by several authors (see, for example, || and references therein). These 'Grishchuk 
gravitons' || are created because the definition of particle states at early and late 
times is not the same and hence an initial vacuum state seen later becomes a multi- 
particle state The direct detection of these waves would provide us with 
a new view of the very early phases of the universe since any graviton should have 
decoupled from matter and from other fields as early as the Planck time (10 _43 s) 
|T|l-|im. In principle, even without detecting them directly, we may obtain useful 
information and constrain important parameters of any cosmological models by an- 
alyzing their indirect effects. In fact, besides probing the current view of the very 
early universe, the presence of a GW background can, for example, cause deviations 
from the Hubble flow for the velocities of galaxies and clusters of galaxies ||, affect 
the timing of orbital motions [|3] and pulsars JTTJ , influence the results of primordial 
nucleosynthesis |2|, f|, and, through the Sachs- Wolfe effect fl^fl , may also contribute 
to the anisotropies of the Cosmic Microwave Background (CMB) JT3], [13 . 



In spite of the importance of the subject, the large majority of the studies on 
graviton production found in the literature deal only with the case of spatially flat 
(e = 0) Friedmann- Robertson- Walker (FRW) universes. In fact, to the best of our 
knowledge, only in two articles the time dependent part of the GW equation is solved 
explicitly for non-spatially flat FRW cosmologies, namely: in the original work of 
Lifshitz ||15|| , where the solution is given for radiation- and matter-dominated universes 
with hyperbolic geometries (e = — 1), and in the recent paper by Allen, Caldwell and 
Koranda |16[ where the authors have considered an elliptical (e = +1) FRW universe 
and a sequence of a de Sitter (inflationary) phase, followed by a radiation-dominated 
era and a matter-dominated (dust) phase. In ref. [|17]] the GW equation for A = and 
any e was established, but not solved for e ^ 0, whereas in the problem of long 
GWs in a homogeneous closed universe (type G3 IX in the Bianchi classification) 
was addressed (see also |19|). 

In the present paper we generalize previous results by studying the creation of 
gravitons in both elliptical (e = +1) and hyperbolic (e = —1) FRW cosmologies. 
For these values of e, we derive the general solution of the GW equation in a model 
where the density p and pressure p of the cosmic fluid obey the equation of state 
^=(7 — l)p, 7<2. We find the solution representing an adiabatic vacuum state 
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2Cfl for a large class of values of 7 (7 7^ | ^n-il > n = 0, ±1, ±2, . . .). By considering 



an arbitrary sequence of cosmic eras with 7 in this class, we evaluate the Bogoliubov 
coefficients which are necessary to obtain the graviton spectrum [El ffl and to give 



an accurate determination of their contribution to the anisotropy of the CMB [|T3 
This exact result implies that it is possible to perform a direct comparison among the 
stochastic background of GW's produced in different FRW geometries (e = 0, ±1), 
and potentially it may also be used to impose constraints on several cosmological 
parameters (the case e = was treated in detail in 0, [|). 

The paper is organized as follows: In Sec. 2 we establish the wave equation for any 
value of e. In Sec. 3 the cosmological model considered is described and the general 
solution for the scale factor in terms of the conformal time, rj, is given. In Sec. 4 we 
derive the general solution for the GW equation in terms of hypergeometric functions, 
whereas in Sec. 5 we obtain the solutions representing an adiabatic vacuum as a 
combination of associated Legendre functions. The mechanism of graviton creation is 
reviewed in Sec. 6 and the Bogoliubov coefficients are evaluated for the model under 
consideration. A brief summary of our results is presented in a conclusion section. 

Unless otherwise stated, the units used are such that h = c = kb = 1. Greek 
indices run from to 3 (spacetime indices) and latin indices from 1 to 3 (spatial 
indices) . 

2 The wave equation 

The study of primordial GW's is done through the use of a perturbation formalism de- 
veloped by Lifshitz ||15|| . This consists in making an expansion in the deviations from 
an unperturbed, non-radiative configuration, and linearizing the Einstein equations 



about that configuration [5T|. The spacetime metric is 

ds 2 = g$ dx^ dx v , (2.1) 

with 

9$ = 9^ + V > ( 2 - 2 ) 

K v < v > (2-3) 

where are small perturbations to the background metric g^ v which can be decom- 
posed into scalar, vector, and transverse-traceless tensor modes. The first two modes 
represent density and rotational perturbations, respectively. The latter transforms as 
a 3-tensor under spatial coordinate transformations, and corresponds to sourceless, 
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weak GW's. We, therefore, consider only perturbations of the Einstein equations 
such that the perturbed value of the energy-momentum tensor, T^, is set equal to 
zero (see Eq. ( |2.20| ) below). 

For a homogeneous and isotropic universe, and in the co-moving coordinate sys- 
tem, the background line element takes the FRW form 

ds 2 = dt 2 - a 2 (t) dl 2 = a 2 (r]) (drj 2 - dl 2 ) , (2.4) 

where t and 7] are, respectively, the cosmic and conformal times, related by 

dt = adr] . (2.5) 

The 3-space metric is 

dl 2 = g ij dx i dx j (2.6) 

= — l — dr 2 + r 2 (d9 2 + sin 2 9 dip 2 ) (2.7) 
1 — er 2 

= d X 2 + F 2 (x) (d6 2 + sin 2 6 dip 2 ) , (2.8) 



Hx) = M^x> = smx e = +1 

= X , e=0 

= sinhx , e = -1 . (2.9) 

We shall use the form of dl 2 in terms of \ , 9 , ip with 0<#<7r, < <p < 2tt, 
< X < 7T for e = +1, and < x < 00 f° r £ — 0, — 1. Note that 

9ij = -a 2 (ty gij . (2.10) 

The Einstein equations read 

R fi v-~Rg fll/ + Ag f tv = -8irGT fXV , (2.11) 

with R^ U1 R, and A being the Ricci tensor, the Ricci scalar, and the cosmological 
constant, respectively. For a perfect fluid the energy- momentum tensor is 

T tl u = (p + p)U^U u -pg^, (2.12) 
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and is the four- velocity of the fluid. The field equations take the form 



S7T 



G Ad 2 



+ 



3 a 2 a 



(2-13) 



^ a a e 

8irGp - A = -2 

a a z a 1 



(2.14) 



where a dot means derivative with respect to the cosmic time t. By assuming that 
the equation of state is 

p=( 7 -l)p , 7 <2, (2.15) 

we get from (|2.13|) and ( (2.141 ) 



P = Po 



a J 



37-2 



ad H 4: — (d 2 + £) — -7Aa 2 = 



1 

2 

As one may check, the above equation has the first integral 



(2.16) 
(2.17) 



<r = { a\Hl + e + ^A a*) (—) ^ - e + I A « 2 



a / 



(2.18) 



where the Hubble parameter is defined by 



H(t) 



a a 



a cr 



(2.19) 



the primes represent derivatives with respect to rj, and p = p(t ), etc. 

To obtain the GW equation we perturb the Einstein equations to first order and 

set 

Sp = 5p = 5U» = , (2.20) 



U» V = • 



(2.21) 



We also impose the gauge conditions 



v„r = o 



(2.22) 



where V v indicates the covariant derivative. Eqs. (|2.21|) and ( |2.22| ) lead to 



h% = 



(2.23) 
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(When a(t) = constant, (|2.22|) constitutes only 3 independent conditions and ( [2.23| ) 
must be imposed separately |T7| .) In the coordinate system of ( |2.4j ) = 5$, and 
Eq. (fE2Ip gives 

V = , (2.24) 

and from (|2~22D 

(2.25) 



Vj h lJ = . 



(Vj represents the covariant derivative with respect to the spatial metric g^.) We 
are then left with only two independent components of W, corresponding to the two 
polarizations of a GW. 



In the gauge used we obtain fT7 | 



d 2 h j dh j If 1 * 



dt 2 



dt a 2 



where 



V 2 V =g lm V t V m W 



J__d_ 

\fg dx l 
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Ira 9hj J ' 

dx m 



(2.26) 



(2.27) 



and g is the determinant of In terms of the conformal time, Eq. Q2.26 ) can be 

recast as 



If we write 



(h/)" + 2-(V')' + lehi = V 2 W . 



x= (x, 0, if) 



(2.28) 



(2.29) 



then it is possible to find mode solutions of fl2.26p (or fl2.28p ), labeled by fc, in the 
form 



hi j (k, t, x) = *(jfe, t) Gi j (k, x) , 
where the functions Gi J (k, x) and ^(k, t) satisfy 

W 2 Gi j (k, x) = -{k 2 - 3e) Gi j (k, x) 
VjGjfcZ) = 0, 
Gi\k,x) = 0, 



(2.30) 

(2.31) 
(2.32) 
(2.33) 



*(M) + 3-*(M) + -e 



(2.34) 
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or still, in the conformal time 



* (k, rj) + 2 — m (k, rj) + (k 2 - e) V(k, n) = 
a 



(2.35) 



In the above equations, k is the co-moving wave number, related to the physical 
wavelength A and frequency uj by 



2ira 

k = — — = uja 
A 



and 15 



< A; < oo , e = 0, — 1 , 
k = 3, 4, ... , e = +1 . 
For e = 0, the symbol stands for 

k = {K- k) , 



(2.36) 

(2.37) 
(2.38) 

(2.39) 



where K represents one of the two possible polarization states of the GW's, i. e., 
K = K a or K = Kf,; k is the co-moving wave vector, 



k = (kijkzjks) , — oo < kj < oo 

1/2 



k = |£| = (l> 4 2 ) 



(2.40) 
(2.41) 



For e = ±1, 



k 
M 
J 



(K; k,J,M), 

-J, -J + l, J, 

2, 3, ife-1 , £ = +1 

2, 3, . . . , £ = -1 . 



(2.42) 
(2.43) 

(2.44) 



(Note that 'our' k 2 coincides with the f3 2 of [0, but is 3£ plus the k 2 of ref. |Tj]], and 
that 'our' is 1 plus the L of ref. ||16|| .) 

It is convenient to express the temporal dependence of the wave in terms of the 
function n(k,t) defined by 

fi(k,t) 



a(t) 



(2.45) 
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Eqs. ( p.34|) and ( |2.35| ) lead, respectively, to 



a 



{k 2 - e) 



H(k,t) = 



(2.46) 



fi (k,rj)+ ( k 2 - e - — J n(k, rj) = . 



(2.47) 



We shall take Eq. ( |2.47| ) as the basic wave equation to be solved. 

The general solution representing the tensor perturbations can therefore be written 

as 

K"(V, = E [Mk) v)G^(k x) + H.C.] , (2.48) 



where H.C represents the Hermitian conjugate of the previous term, and J2 refers to 



summation over discrete values and integration over continuous ones, i. e. |17L 22 



E 

k 



E 

K 

E 

k,J,M 



E 

J,M 



/roo 
d 3 k = J2 k 2 dk , e = 

is. K 

E E 



K k,J,M 



E L dk E 

J,M 



K 



K b 

E , 

K=K a 
oo k—1 J 

EE E 

fc=3 J=2 M=-J 
oo J 

E E 

J=2 M=-J 



e= +1 



e = +1 



-1 . 



(2.49) 
(2.50) 
(2.51) 
(2.52) 



In the coordinates used, the spatial functions are specified by Eqs. ( |2.31|) -( p.33|) 
and by 

<V = . (2.53) 
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We shall impose the normalization condition 

G/(jfe, x) G*^(k!,x) sfgtfx = 16 tt G 8(k, k') 



(2.54) 



where the asterisk indicates the complex conjugate and 5(k, k') is the delta function 
with respect to HI 



1. e., 



^g(k')6(k,k')=g(k) 



(2.55) 



It proves convenient to choose the phases of the G^ u (k, x) in such a way that 

G;»(k,x)=G^(-k,x), (2.56) 

where 



k = -(K; k) = {K;-k) , e = , 

k = —(K; k, J, M) = (K; k, J, —M) , e = ±1 . 



(2.57) 
(2.58) 



(See [f22[] for details.) 

In the classical theory A(k) and A^(k) in Eq. ( 2.48|) are just complex constants 
(Fourier-Bessel coefficients), whereas in the quantized theory they become annihila- 
tion and creation operators, respectively. As usual, these operators are required to 
satisfy the commutation relations 



[A(k),A*(k')] = 5(k,k') 
[A(k),A(k')} = 0. 



(2.59) 
(2.60) 



We further impose the following normalization condition on the Wronskian of the 
time-dependent part of the unified solution 



or, equivalent ly, 



^(k,ri)^*'(k,rj) - V*(k, rj) rf) 



H(k,r)) n*'(k,T)) - fi*{k,r}) f/(k, 77) = i . 



(2.61) 
(2.62) 
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3 The solution of the field equations 



3.1 General equations 



The problem of determining the spectrum of primordial GW's in a consistent way 
may be summarized in four basic steps: 

1. First of all one finds the a(rf) for the cosmological model under consideration 
by solving ( p. 17] ) or, equivalently, Q2.18|) . This is required in order to obtain 
the explicit form of the wave equation (|2.47|) . Particular solutions of the field 
equations have usually been considered in the literature. In this work a general 
treatment will be presented. 



2. The solution of ( |2.47| ) should then be found. As remarked earlier, some approx- 



imate methods have been used in the literature (see, for example, f2~3|). 



3. By studying the behavior of the solutions in the high-frequency limit, the partic- 
ular solution of (|2.47|) representing an adiabatic vacuum state |2(| is identified. 



4. Finally, the Bogoliubov coefficients are evaluated using, for instance, the general 
procedure described in [fjj. 

In this section we deal with the first step by writing the solutions in a convenient 
way for our purposes. Following the procedure developed in ref. pl| , we initially 
recast ( 2.17 ) using the conformal time 



a (37 — 4) (a 
a 2 \ a 



i\ 2 



37-2 1 A 2 

+ — -e 7 A a 2 

2 2 ' 



(3.r 



For 7 = 2/3 we make the substitution 



b = In a 



(3.2) 



to get 



which can be easily solved for A = 0. The scale factor is then 



(3.3) 



a( v ) =a e aoHo{v - Vo) , 7 = - , A = 



(3.4) 
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For 7 ^ 2/3, by replacing 

b = a 1/q , (3.5) 



« s »H' (3 - 6) 



Eq. (|3.1|) takes the following form 



b" + -h-^b 2 « +1 = 0, (3.7) 



which formally is the nonlinear equation satisfied by a forced harmonic oscillator. 

From now on, we restrict our analysis to models with vanishing cosmological 
constant. In this case, ( p.7| ) reduces to the equation of a simple harmonic oscillator, 
whose solution is 

b=c l em(J^ri + c^ , 7 ^jj , A = 0, (3.8) 

and from ( |3.5| ), we write the unified expression for the scale factor as 

/ sin 6 \ 9 , 2 , . 

<V ) = ao K^T ' ' A = ' 3 ' 9 

VsmBo/ 3 



where, 



e = ^(77-770) + 8 , (3.10) 
m 

60 = 6(r/ ) = arctan ij- j , (3.11) 



and by definition 

ao = a(r) ) , (3.12) 

Ho ee #(r/ ). (3.13) 

Note also that the unified expression for the Hubble parameter is given by 

H( V ) = ^^ £COte . (3.14) 
\q\ a 
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/From flnSD with A = and 
the density parameter is then 



3H 2 



(3.15) 




1 + 



e 1 



(3.16) 



a- 



2 H 2 (cos6) 2 ' 



In the limit e — » 0, ( [3.9|) reduces to the general form of the solution for the spatially 
flat case (Eq. (28) of ref. H). Since the graviton spectrum for e = was studied in 



3.2 Elliptical models (e = +1) 

In these models, B and ©o are real and, for q > (7 > 2/3, non- inflationary sce- 
narios), Go = arctan(ao-ffo), whereas for g < (7 < 2/3, inflationary scenarios), 
60 = 7r — arctan(ao-f^o)- Therefore, for any q, < Q < n and sin 60 > 0. For a(rj) 
to be positive we must then require < < n. 

For q > (a < 0) there is an initial singularity at = (rj s = r/ — g©o). 
An expansion phase occurs for < < 7r/2 until a(rj) reaches the maximum value 
a- m ax = ^/(sinOo) 9 at = 7r/2 (7] max = rjo + g(f-0 o J). This is followed by 
a contraction stage for (tt/2) < < n, and finally by a 'big-crunch' at = 7r 

{fjs =V0 + <?(?!" - ©o))- 

For q < (a > 0), a is arbitrarily large as — > (77 — » Tjooi = Vo — \q\®o), 
decreases for < < 7r/2, reaches the minimum a min = a (sin ©o)' 9 ' at = 7r/2 
(^mm = Vo + M (f — ©0)), grows for (7r/2) < < 7r, and becomes arbitrarily large 
as -> 7r (r? -> 77003 = ?7o + |g|(7r - O )). 



Note that cos© is positive for (0 < © < f , q > 0J, f| < © < n , g < 0J , and 
negative otherwise. 

3.3 Hyperbolical models (e = —1) 

For e = — 1 we define 



detail in @, §], in what follows we shall treat only the cases e = ±1. 








1 



(3.17) 
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$0 = = arcth (aoH ) , 

q 1 



(3.18) 
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where |ao-f^o| > 1- Hence, 



sinh $ 



E = t- 



sinh (M $0 ). 
q coth $ 



a 



\q\ sinh$ 
q sinh $ 



191 a 



(cosh*)' 



(3.19) 

(3.20) 
(3.21) 



For an expanding universe we must have 

a H > 1 . (3.22) 

Consequently, $ > 0, sinh$ > 0. The condition a > then requires $ > 
(V > Vo — Q ) f° r Q > 0> an d $ < (77 < r] + |g| $ ) f° r Q < 0- 

For g > there is a singularity at $ = (r/ s = rjo — q $o) an d a — > oo as $ — > oo 
(77 — > oo). For g<0, a-»0as$-» — oo (rj — > — oo) and a — > oo as $ — > 
(77 -> r/o - g$ )- 



4 General solution for the wave equation and the 
effective potential 

For the cosmological model described in Sec. |3|, the wave equation ( |2.47| ) reads 



jj, (k,rj) + (k - e - a H ) fi(k, rj) = 



,7 3 , 



(4.1) 



A; 2 + 



g (sin©) 2 



Kk,v) = » 7 ^ 3 ' 



(4.2) 



The solution of ( |4.1| ) is trivial and we shall concentrate only in the case 7 7^ ~. 

Note that, as it happens in the spatially flat geometry, the effective potential 

' ;i-*)e 



q (sin 0) 2 



(4.3) 
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approaches zero as 7 — > 4/3. If we define the Hubble co-moving wave number by 



2 7r a 2tx a 



A 



H 



H 



and 



S = 



eff 



k 2 ' 



we see that the ratio S is time dependent, i. e., 



(4.4) 



(4.5) 



S(e(r,)) 



An' 



1-q 



(cos6) s 



(4.6) 



Therefore, the amplification condition |2|, |], |2B 

k 2 < y e// 



(4.7) 



does not necessarily represent modes outside the Hubble radius, A ^> A# (k 2 kjj). 
In fact, this will happen only for those time intervals where 5(0(77)) < 1- I 11 elliptical 
geometries, however, we may have S(Q(rj)) > 1 as long as 



1 



1-q 



q 



< cos < 



2tt \ 



1-q 



q 



(4.8) 



that is, close to a max for q > and close to a m i n for q < 0. As | cos0| < 1, ( |4.8| ) is 
restricted to hold for q in one of the following intervals: (1 + 47r 2 )~ 1 < q < 1, q > 1, 
or g < (1 - 47T 2 )- 1 . 

In the hyperbolic geometries, 5(0(r/)) > 1 if 



cosh<3> < 



2tt \ 



1-9 



(4.9) 



i. e., near the singularity when q > 0, or as a — > 00 for q < 0. This can happen only 
for < q < (1 + 47r 2 )~ 1 or (1 - 4 7r 2 ) -1 < q < (as cosh$ > 1). 
In order to solve (14.21) we first make the substitution 



(sin9)' 



(4.10) 
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and write 



MM) =(7) f( z )- 



(4.11) 



It is easily seen that the function /x defined by ([4.11 ) is a solution of the wave equation 
as long as / satisfies 



d 2 f 



q + 



(4.12) 



which is a hypergeometric equation |27| with parameters 



B 



C 



? 1 



(4.13) 
(4.14) 
(4.15) 



The general solution of the gravitational wave equation is then given by Q4.ll ) 



with f(Z) being the general solution of the hypergeometric equation ( 4.12 ). The 
explicit form of this solution depends on several relations involving the parameters 
A, B, C, and may be found in p^-pfl. 



5 Adiabatic vacuum solutions 



The next step is to find the particular solution which represents a 'physical' vacuum 
state. As it is usually done in the literature, we adopt the adiabatic approach for 
defining particle states. In the high-frequency limit the adiabatic definition reduces 
to the usual positive-frequency Minkowski modes (~ e~ lkv ). A detailed discussion of 
this prescription can be found in |^0| and, in our model, it can be easily accomplished 
whenever 

C^O, ±1, ±2, ... , (5.1) 

which is equivalent to 



, 1 

Q T 1 n 

HT 2 



, n = 0, ±1, ±2, ... , 



(5.2) 
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or still 



, 2 (2n + 1) 



3 (2n - 1) 



n = 0, ±1, ±2, ... . 



(5.3) 



Under this restriction (which covers interesting cases such as 7 = 0, 7 = 4/3, 7 = 1), 
the general solution of Eq. ( |4.12| ) can be written as ||27||— 



f(Z) = d F{A, B, C- Z) + C 2 Z 1 ' F{A', B',C';Z), (5.4) 
where F stands for hypergeometric functions, C\ and C 2 are arbitrary constants, and 



A 

B' 

d 



A-C + l, 
B-C + l, 
2-C . 



(5.5) 
(5.6) 
(5.7) 



The solution for the wave equation can then be expressed as 



z y/2 



fi(k, v)=C ijj [F(A, B, C; Z) + Co Z l - C F(A , & , C ; Z) 



4 ' is-*' 



(5.8) 



The constant Co must be found by imposing the adiabatic constraints, whereas C is 
determined up to an irrelevant phase by the normalization condition fl2.62| ). 
The fact that for both values of e, 



C = A + B 
C' = A' -B ■ 



l 



2 ' 
l 



(5.9) 
(5.10) 



enables us to write the hypergeometric functions appearing in ( |5.8|) in terms of asso- 



ciated Legendre functions |2"T||-f2"9[ . 

We anticipate the final result derived below and write the adiabatic solution for 
the GW equation in the unified form 



fi(k,r}) 



T(m + l) 



-m) 



11/2 



4/r(m + l) 



(sin6) 1/2 £(X) 



(5.11) 
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where: 



c(x) = v- m {x) _ e T^(^)£0_^ OT) (512) 

1 [y + m + 1) 

X = |cos6| , (5.13) 

1 3 / 2 -7 \ , 
mS «-2 = 2^j' (5 ' 14> 

v = y/e\q\k--, (5.15) 

/ - Im(e-^m ^~ m + 1 H , (5.16) 
\ T(z/ + m + 1) J ' V ; 

Im{z} stands for the imaginary part of z, ip is a constant phase (whose value is 
irrelevant for the evaluation of the Bogoliubov coefficients), T is the gamma function, 
and V™ represents the associated Legendre function of the first kind, of order m and 
degree v. The upper (lower) sign applies if cos 6 > (< 0). 

Note that for e = —1, cosG = cosh$ > 1, but for e = +1, cos0 > during an 
expansion phase if q > and during a contraction phase if q < 0. We also have 

0<Z<1,0<X<1 (£ = +1) (5.17) 

and 

- oo < Z = -(sinh $) 2 < , 1 < X = cosh $ < oo (e = -1) . (5.18) 



In ref. ^ZE\ the 'slanted' symbols P™(X), Q™(X) are used to indicate the asso- 
ciated Legendre functions of the first and second kind defined for Re{X} > 1. The 
'straight' symbols P™(X), Q™(X) are used for \X\ < 1. We adhere to this notation 
whenever we are dealing explicitly with the cases e — — 1 and e — +1, respectively. 
The calligraphic symbols V™ and will be used only in unified expressions like 



( |5.11| ) that are valid for both values of e. (Note that, as it was remarked in \TB\, there 



is a misprint at page 999 of ref. |28|: the first paragraph should refer to 'straight' P 
and the second paragraph to 'slanted' P.) 

Some functional relations involving the associated Legendre functions (relations 
with the hypergeometric functions, asymptotic behavior, derivatives, etc.), which are 
needed to derive the adiabatic solution and the Bogoliubov coefficients, are distinct 
for P™(X) and P™(X) p8[. We must therefore treat the two cases separately. 
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For the elliptical models, and due to equations ( |5.9| ) and (|5.10|) , we can use Eq. 



(15.4.13) of ref. to write F{A, B, C; Z) and F(A' , B' , C'; Z) in terms of P~ m (X) 
and P™(X), respectively. If we require /i(fc, rj) to have the desired asymptotic behavior 
for large k, the constant Co may be determined. This can be done through the use 
of Eq. (8.10.7) of [p7 |. In order to evaluate C we apply the condition ( |2.62| ) on the 



Wronskian of the solutions by using Eq. (8.741.1) of p8| . Whenever we are forced to 
deal with P™(— cos 6), we use Eq. (8.737.2) of [28]. The final result is then given by 
(PH) - (P^) with e = +1 and V™ i-> P™. 



For the hyperbolical geometries, equations ( |5.9D and ( |5.10|) enable us to use Eq. 
(15.4.12) of [|7| to write F(A, B, C; Z) and F(A' , B',C'; Z) in terms of P~ m (cosh $) 
and P™(cosh <£>), respectively. To analyze the asymptotic behavior for large k we make 
sucessive use of equations (8.723.1) of |2pJ and (15.7.1) of [0. The constant Co is then 
determined. For the Wronskian condition, we first write P~ m (X) as a combination 
of F?{X) and Q™{X) using Eq. (8.736.1) of @ and then apply (8.1.8) of |Z7|. We 



finally get (g^)-fUg) with e = -1 and ^ P™. 

6 The Bogoliubov coefficients and the graviton spec- 
trum 

6.1 Review of the graviton creation mechanism 

The general procedure to calculate the graviton spectrum was explained in detail in 
0. We shall give here only a brief resume of this method and then evaluate the 
Bogoliubov coefficients for our present model. 

In realistic cosmological models the equation of state (|2.15|) will change its form 



during the cosmic evolution. Prior to a time r\ r the equation reads 

P(r-l)(v) = (7r-l - 1) P(r-l){v) , (V < Vr) , (6.1) 

whereas for 7] > r] r 

P(V){V) = (7r - l)P(r)fa) , {V > Vr) ■ (6.2) 

We are working in the sudden transition approximation where the change in 7 occurs 
instantaneously at the transition time t] r . (See for the conditions of applicability 
of this approximation and other details.) It is possible to impose the continuity 
conditions, 

a (r _i)(?7 r ) = a(r)(r) r ) , (6.3) 
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a'(r-l)(f7r) = a[ r )(Vr) ■ (6.4) 

(9(r)(v) indicates the appropriate form for the quantity g in the interval [rj r , r) r+ x\. 
The transition times are T] r , t] r+ i, etc. For the first phase we take rjo to be any instant 
in that phase. The notation g r indicates g( r )(Vr), etc.) 
We write 



y(r-i)(k,r)) 



1 



a (r) (rj) 



a r (k)p(r)(k,rj) + p r (k)fi{ r) (k,r)) , r? > r] r , (6.5) 



and 



G(r)(^) 



a£(fc) /i( r _i)(/c, 77) - V* r -i)(k, ^) 



, 77 < ?7 r 



?7 > ?7r 



(6.6) 



where fi( r -i) and /i( r ) are solutions to the wave equation describing an adiabatic 
vacuum in their respective regions. The two modes !V( r -i) and ^(r) define two quan- 
tizations of the field associated with two different Fock spaces, but the associated 
operators in each case represent physical particle observables only inside their respec- 
tive eras. It is easily seen that 



3> (r _i)(fc,77) = a r (k)y (r) (k, V ) + (3 r (k)yt r) (k,ri) . 
Due to ( |2.62| ), the Bogoliubov coefficients [p0| , a r (k), /3 r (k) satisfy 



(6.7) 



\a r (k)\ 2 - |ft(fc)| s 



(6.8) 



and the operators in the two epochs are related by a Bogoliubov transformation 



A {r) {k) = a r {k) A (r _D(fc) + P* r (k) A\ r _ l} {-k) 
V-i)(*) = a* r (k) A (r) (~k) - #(*) A\ r) (-k) . 



(6.9) 
(6.10) 



Consequently, the vacuum state at the region (r — 1), labelled |0( r _i) >, is annihilated 
by Ai r -i) but not by Am, and, if N^(k) is the number operator for the mode k at 
stage (r), 

AT (r) (£)=A| r) (£)A (r) (£), (6.11) 
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then [g [20 1 

< N{k) > r =< (r _ 1) |iV (r) (fc)|0 (r „ 1) >= \p r (k)\ 2 . (6.12) 

Hence the vacuum state |0( r _i) > is a multi-particle state when we use the definition 
of particles appropriate for the epoch rj > rj r . This is interpreted by saying that 
particles have been created by the expansion dynamics. 

In a multi-stage model we need to relate the operators A and of stages separated 
by several transitions. This can be done recursively by H, 0L 19J 



«T r (fc) = a r (k) ar^^k) + (3*(k) (3 Tr ^{k) , 



(6.13) 
(6.14) 



The Bogoliubov coefficients can be found if we impose the continuity of 3^(r-i) and 
y^r) and their first derivatives at r\ r . The result is 0, f| 



a r {k) 

/3 r {k) 



fi* r) (k, rj r ) /i' (r _ 1} (k, r] r ) - /z (r _i) (k, rj r ) nf r) (k, rj r 

M(r-l) (k, T}r) /i' (r) (k, rjr) ~ fifr) (fc, T) r ) M( r _l) (fc, 



(6.15) 
(6.16) 



The power spectrum of the graviton background can be described by the quantity 
P g (u), defined in such way that P g {uj)duj represents the energy per unit volume 
between the frequencies u and u + du [0, f|, [|. In the units used in this paper 



Pa(w) = ^ < N L 



7T 



w(fc) 



> 



(6.17) 



The above expression is valid for waves with X < Xh- Modes outside the Hubble 
radius (A > Xh = H~ r ) cannot contribute to the energy density since this is a 
locally-defined quantity (see @, |, g, g, JZQ for details). 

Note that equations (|6.15|) - (|6.17|) are independent of the cosmological model con- 
sidered. Eq. Q6.12 ) assumes that no particles were initially present. If the initial state 
is not the vacuum, then 'stimulated' graviton creation may also occur 
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6.2 The Bogoliubov coefficients 



In our model, the expressions for the Bogoliubov coefficients evaluated from ( |5.11| ), 
( |6.15| ), ( |6.16| ) may be simplified if we use the continuity conditions ( |6.3| ), (|6.4|) , which 
read 



sin 9 



(r-l)(Vr 



sin 6 r _i 



1r-l 



a r -i 



(6.18) 
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cote (r _i)(r7 r ) 



q r -i \q r \ 

For e = —1, (|6.18|) and ( |6 . 1 9|) reduce to 



|g r _l| _5V_ CQt q = |gr-l| «r #r 



5r-l V £ 



sinh$ (r _ 1) (?7 r ) 

coth$ (r _ 1) (?7 r ) 
For e = +1 the final result is 



gr-1 
|?r-l| 



sinh $ r _i 



?r-l 



— coth $ r = — r —^— a r H, 



a r _i 

|?r- 



a r (fc) 



i e 



i (lp r -!-1p r ) 



q r -i 



i e 



i (4> r -l+4> r ) 



q r 

q r -i 



q r 



1/2 



1/2 



q r -i 



Br-iB r . (Y r \ 1 ' 2 



cos r | = sin r |a r if r | = Vl r , 



sin 9 r _i / a r \ 
cosB( r _i)(r7 r )| = — : — - — I I X r 



sinB r \a r -i 

|r(— m r )\ T(u r — m r + 1) 
r(— m r ) T(v r + m r + 1) 



sin m, 7r 



1 V2 



2 2 |g r _ x | |g r | 



2 2 |g r _ x | |g r 



Q£) +1 (r r ) P™;(x r ) - Q^(Y r ) Pt r)+ \x r ) , 



?r-i 



9r-l 
9r-l 
9r-l 



p<U+i(y r ) P ™ ;( x r ) - P^) P^ )+1 (X r ) . 
P^J +1 (Y r ) Q^(X r ) - P l ;(Y r ) Q^ +1 (X r ) 



(6.19) 



(6.20) 
(6.21) 



(6.22) 

(6.23) 
(6.24) 
(6.25) 

(6.26) 
(6.27) 
(6.28) 
(6.29) 
(6.30) 
(6.31) 
(6.32) 
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l r = m r _i , (6.33) 
n r = z/ r _i . (6.34) 



The quantities 6 r , ip r , q r , m r , u r , I r are those defined by equations ( |5.11 



( |5.16j ) appropriate to the epoch r] r < rj < rj r+ i (and similarly for r _i, etc.). The 



upper (lower) sign applies if cosO r > (< 0). 
For e = — 1 the coefficients can be written as 

a r (k) = ie ^-i>r) e -i{ir+mr)-K ( ^zl)^ ^zl B ' sinh $ r (lA ' JJ a ,(6.35) 

V Qr J n \X r J 

fn \ 1/2 R R fV \ 1/2 

(3 r (k) = ie (^- 1 +A) e - l (ir+m r ), IQr-i\ ^zI^H sinh $ r ±M ^,(6.36) 

V Qr J n \X r J 

X r = cosh$ r , (6.37) 

Y r s c sh* {r _ 1)( „) = (^) 91 X r , (6.38) 

D ar = + E r Df ] , (6.39) 

^ = , (6-40) 

gim.Tr i m {r*(z/ r -m r + l)r(z/ r + m r + l)} 

£ r = in— — — ■ — - — , (6.41) 

sm(m r 7TJ 1 (z/ r . — m r + 1) 1 \y r + tHt + 1) 

E? } = ^Q£ )+1 (^)Qr(X)-Qt(nOQ(r )+1 (X r ), (6.42) 

<Zr-l 

Tjf = J^Q(M+i(F r )P-(X r )-Qt(F r )P^ )+1 (X). (6.43) 



7 Concluding remarks 

We have addressed the problem of graviton creation in elliptical (e = +1) and hy- 
perbolic (e = —1) Friedmann- Robertson- Walker cosmologies with a gamma-law type 
equation of state. After establishing the general solution for the scale factor a(rj) 
we have obtained the general solution for the gravitational wave equation. We have 
then found those modes representing an adiabatic vacuum state for 7 7^ | ff^zj), 
n = 0, ±1, ±2, . . .. This includes most cases of physical interest, such as a de Sitter 
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phase (7 = 0), a radiation-dominated era (7 = 4/3), and a dust phase (7 = 1). By 
considering an arbitrary sequence of epochs with 7 in this class, we have derived the 
Bogoliubov coefficients associated with graviton creation in this model. These coeffi- 
cients are necessary to obtain the gravitational wave power spectrum, P g (u), and to 
give an accurate determination of the graviton contribution to the anisotropy of the 
CMB 0. 

The application of our results to a definite model, that is, to a definite sequence of 
cosmic eras corresponding to particular values of 7, can, in principle, be used to com- 
pare the behavior of the relic gravitons in the three possible geometries (e = 0, ±1). 
In fact, by imposing the observational constraints associated with the anisotropy of 
the CMB, with the nucleosynthesis of 4 He, and with the limits on the present value of 
Q, it may be possible to obtain theoretical constraints involving several cosmological 
parameters, such as the value of 7 in an initial inflationary period, the value of H at 
the end of inflation, and Q itself. 

Some other interesting questions also arise. For example, in elliptical models, 
what happens to graviton creation as the expansion maximum is reached? And, for 
models with a big-bounce (q < 0), what can be said about this production near the 
expansion minimum? |33| Some of these issues will be addressed in a forthcoming 



communication [j3l| . 
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